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abstract 


This  report  discusses  properties  and  use  of  a  highly  resolving 
r.f.  cavity  probe  system  for  wake  trails  in  hypervelocity  model  ranges. 
Data  are  taken  from  two  distinct  modes  of  the  system  and  yield  electric 
susceptance  at  an  average  axial  position  inside  the  cavity.  This  is 
important  in  wake  trail  study,  since  it  ties  together  a  wake  electrical 
measurement  with  distance  behind  a  hypervelocity  model. 
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NOMENCLATURE 


aiJ 

b 

d  =  e  +  f 
IE,  E 
h 

ft,  H 

k 

K-l 


unsubscripted,  wake  radius 
subscripted,  matrix  element 
unsubscripted,  cavity  radius 
subscripted,  matrix  element 

cavity  length  in  terms  of  rearward  e  and  forward  f  division, 
electric  field  vector  and  scalar  component 
propagation  constant 

magnetic  field  vector  and  scalar  component 

wavenumber 

susceptibility  in  e 


K-l 

M-l 

q 

Q 

r,* 

S  =0+  JCD 


TM  ,TM  . 
mn  mnl 


averaged  susceptibility  in  d  =  e  +  f 
susceptibility  in  f 
normalized  Q  factor 
circuit  Q  factor 

radial  coordinate  and  unit  vector 
complex  radian  frequency  u>  and  decay  rate 
respectively  waveguide  and  cavity  TM  waves 
variable 


y  =  g  +  j(u£  -  ~)  normalized  admittance  in  Section  II 


y 

z  ,£ 


variable  in  Appendix  A 

longitudinal  coordinate  and  unit  vector 
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7 


TIT,  n 

a 

<J> 

* 

¥ 

CJ 


propagation  constant 
unsubscripted ,  increment 

subscripted,  inner  product  of  unit  Bessel  functions  Ms 
eigenvalue 

vector  potential  and  scalar  component 

time  decay  rate 

radial  component  of  potential 


Ml 


=  1 


unit  Bessel  function, 
longitudinal  component  of  potential 
radian  frequency 
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SECTION  I 

INTRODUCTORY  REMARKS 


This  final  report  discusses  properties  and  use  of  a  highly 
resolving  r.f.  cavity  probe  system  for  hypervelocity  wake  trails. 

Data  are  taken  from  two  distinct  modes  of  the  system  and  yield 
electric  susceptance  at  an  average  axial  position  inside  the  cavity. 
These  are  important  data  in  wake  trail  study,  since  they  tie  wake 
electrical  measurement  together  with  distance  behind  a  hypervelocity 
model.  Below,  following  brief  description  of  a  current  trend  in 
wake  trail  study,  the  basic  probe  system  will  be  indicated,  and  an 
outline  made  of  subsequent  consideration  of  design  problems  of  the 
multimoded  probe. 

Shock  dominated,  wakes  of  blunt  bodies  in  air  are  characterized 
by  rapid  decay  of  large  initial  ionization.  Decay  is  more  rapid  as 
ambient  pressure  increases,  and  ionization  is  quenched  by  the  nearly 
complete  reaction  of  ablated  material  with  air  --  when  ablation  is 
present  —  and  reassociation  of  air  components.  A  rough  picture  of 
the  reacting  sequence  (following  Kornegay ,  Ref.  8)  is  given  at  least 
in  the  following  steps:  electron  dissociation  in  air,  electron 
association  with  oxygen,  completed  by  ionized  oxygen  attachment  to 
positive  ions.  The  data  of  Kornegay  are  similar,  except  for  ioniza¬ 
tion  level,  to  the  data  of  Hendrix  and  Mayfield. 

The  process  can  be  seeded,  most  likely  by  ablated  material  or 
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by  contaminants  of  air.  It  is  known,  for  example,  that  traces  of  sodium 
in  air  can  contribute  significantly  to  electron  population  and  to  collision 

i 

frequency  because  of  extremely  low  ionization  potential  (Refs.  2  and  6). 
Such  behaviour  of  contaminants  provides  sources  of  electrons  to  speed 
the  pure  air  reactions.  However,  it  appears  that  data  on  initiation 
of  ionization  for  these  chemical  processes,  and  on  the  role  of  air 
contaminants  therein,  have  not  much  been  measured. 

In  all,  air  is  a  complicated  substance.  Its  properties  are  often 
better  measured  experimentally  than  predicted  theoretically.  It  is  the 
purpose  of  this  report  to  describe  characteristics  of  a  wake  trail 
measuring  system  (a  multiple  moded  r.f.  cavity)  which  retains  two  useful 
properties  in  its  measurements:  namely,  continuousness  and  high  resolution 
of  data,  axially  along  the  wake  trail.  In  a  sense,  the  system  is  a 
prototype.  It  could  also  yield  resolved  radial  data. 

A  primitive  block  diagram  for  the  system  is  shown  in  Figure  1. 

There  are  other  approaches  to  the  measurement,  of  which  this  is  the 
simplest  one.  For  definiteness,  it  is  assumed  excited  in  the  TMq2q 
mode  with  output  by  the  TMq^  mode.  After  calibrating  without  the 
wake,  only  two  readings  are  needed  when  the  wake  is  present,  the 

i 

reflection  and  transmission  coefficients  or  equivalents.  The  amplifier 
inserted  in  the  output  side  will  help  accuracy  of  the  measurement. 

Accuracy  is  best  when  input  and  output  both  are  read  at  about  the  same 
level.  The  degree  of  accuracy  in  the  design  calculations  done  here 
can  be  illustrated  by  an  example.  If  wake  Q  is  about  100  (exaggerated). 
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then  the  system  will  yield  a  complex  dielectric  constant,  averaged 
only  radially,  at  a  point  inside  a  0.1  inch  interval  of  a  10  inch 
long  cavity. 


Fig.  1  System  Block  Diagram 


It  is  well  to  point  out  that  multimoding  the  system  requires  only 

that  reading  be  done  on  two  different  modes .  It  is  not  necessary  that 

these  be  simultaneously  resonant.  It  is  merely  convenient.  It  can 

retain  small  size  of  equipment  when  modes  resonate  together.  There  is 

however  a  maximum  of  three  modes  that  can  resonate  at  once  in  the  cavity. 

Thus  a  higher  moded  system  would  use  separate  cavities  or  separate 

frequencies  to  operate.  The  choice  can  be  one  of  practical  convenience. 

What  is  necessary  for  system  operation,  however,  is  just  measurement 
* 

on  distinct  modes. 


*Even  this  is  not  unqualified.  The  same  modes  in  different  size  cavities 
could  also  yield  a  reading  of  wake  susceptance. 
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Material  discussed  subsequently  is,  in  outline,  as  follows. 
Section  II  will  present  the  system  equivalent  circuits  and  consider 

i 

the  relation  of  circuit  resonances  with  measured  circuit  parameters. 
Results  are  shown  in  Table  II.  In  Section  III,  by  solving  electro¬ 
magnetic  field  equations,  circuit  resonances  and  measured  parameters 

i 

are  related  to  cavity  and  wake  physical  properties.  These  results 
are  shown  in  Table  VII.  The  interesting  question  of  optimizing  the 
wake  measurement  is  considered  in  Section  IV.  One  is  asked,  equiva¬ 
lently,  to  measure  the  wake  inside  the  cavity  at  an  axial  position 
that  divides  the  weight  of  ions  in  half.  These  results,  in  Table 
VIII,  locate  the  wake  average  susceptance  at  its  average  position 
inside  the  cavity.  An  appendix  discusses  error  (the  largest  one) 
in  solving  the  field  equations  in  Section  III. 
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SECTION  II 

PERTURBED  CAVITY  EQUIVALENT  CIRCUIT  PARAMETERS 

The  cavity  can  be  viewed  as  a  network  nested  inside  another 
network  called  the  internal  circuit.  This  is  shown  in  Figure  2. 

Whatever  attaches  to  the  ports  of  the  internal  circuit,  one  can  call 
the  external  circuit.  An  ordinary  analysis  would  determine  proper¬ 
ties  of  the  cavity  by  solving  an  electromagnetic  boundary  value 
problem,  but  neglecting  wall  losses.  Wall  losses  are  then  usually 
either  easily  measured  or  else  easily  calculated  from  the  now  known 
lossless  wall  solution.  All  this  would  determine  the  elements  of 
the  internal  circuit  of  Figure  2.  The  unit  conductances  represent 
cavity  wall  losses. 

In  case  the  cavity  is  loaded  with  a  slightly  lossy  filling, 
then  the  lossless  solution  can  be  made  to  serve  here  also.  This  is 
true  even  when  variational  approaches  are  taken,  so  long  as  system  Q 
factors  are  high,  as  they  are  assumed  to  be  for  the  present  cavity. 

All  the  above  mentioned  procedural  step$}  will  be  taken  in  considering 
this  design  of  a  two  moded  cavity.  However,  the  present  section  will 
focus  on  behaviour  of  the  circuit  needed  to  represent  the  system. 

An  extended  analysis  would  account  also  for  the  coupling  circuits 
to  the  cavity,  but  these  additional  loadings  are  ineffective  for  the 
most  part,  and  are  overlooked  therefore.  They  could  be  considered 
separately.  Sources  to  drive  the  cavity  and  meters  to  measure  transmissions 
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Fig.  2  Internal  and  External  Circuits 
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and  reflections  are  located  in  the  external  circuit.  In  order  that 
sources  and  metering  equipment  not  contribute  reflections  to  the  system, 
beyond  those  due  only  to  the  internal  circuit,  external  circuits  are 
assumed  adjustable  to  match  with  g  =  1. 

The  cavity  circuit  exists  in  two  states,  measurements  on  both 
of  which  are  required  in  order  to  characterize  the  wake.  The  state 
shown  in  Figure  3(a)  occurs  for  the  empty  cavity.  In  this  state,  the 
circuit  parameters  Lq,  Cq,  the  double  resonance  measurement  frequency 


and  the  numbers  qQ  and  q  can  all  be  found  or  measured.  Here, 


q 


o 


2i 

*1 


—  (2) 


is  the  ratio  of  internal  circuit  Q's  (TMqj^  to  TM^q)  of  the  empty 
waveguide,  while 


q 


=  q. 


—  (2a) 


where  and  k2  are  wave  numbers  respectively  for  TMgg  and  TMq^  circular 
waveguide  modes. 

In  its  second  state,  that  shown  in  Figure  3(b),  the  cavity  holds 
a  model  wake.  It  is  here  that  results  of  an  electromagnetic  boundary 
value  problem  will  be  used  to  relate  to  circuit  parameters.  However,  the 
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wake  is  present  for  only  a  few  hundred  microseconds  and  this  means  there 
is  no  time  to  check  out  completely  the  Figure  3(b)  circuit  (against 
frequency).  Instead,  a  quick  measurement  is  made  on  the  circuit  in 
Figure  3(c) .  By  inference,  this  is  related  to  the  Figure  3(b)  circuit 
from  results  of  the  solved  electromagnetic  field  equations.  The  Figure 
3(a)  circuit,  the  first  state,  having  been  checked  out  before  the  model 
flight,  then  all  parameters  in  the  equivalent  cavity  circuit  therefore 
would  have  become  measured . 

The  cavity  circuit  admittance,  when  the  wake  is  present,  has  two 
zero  resonance  conditions, 

y-  =  0  —(3) 

and 

qyj^  +  (l  +  q)y2  =  0  »  —  (4> 

whether  looking  in  at  Side  1  (TMq20  terminals)  or  at  Side  2  (TMq-q  terminals). 

This  reflects  the  property  of  the  wake  that  it  does  not  merely  shift  the 

individual  modal  frequencies,  but  rather  it  couples  them  together  making 

resonances  such  that  the  modes  are  jointly  present.  The  fact  that  the 

same  split  and  shift  in  resonance  occurs  at  Side  1  as  at  Side  2  follows 

* 

from  solution  of  the  field  equations. 

It  is  the  same  effect  that  causes  the  change  of  q  to  q  between  the 
Figure  3(a)  and  the  3(b)  and  3(c)  circuits.  It  is°due  to  the  fact  that 
cavity  magnetizing  current  transfers  from  terminals  1  to  2  by  the  real 
number  factor  jk^l2  different  from  the  transfer  from  terminals  2  back  to  1. 

Ivl 
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By  contrast,  since  the  cavity  circuit  of  Figure  3(b)  is  not 
symmetrical,  the  admittance  pole  resonances  are  not  constrained  to  be 
the  same  at  both  terminals.  Here,  the  pole  resonance  condition  on 
Side  1  is  that 

qy{  +  y2  =  o  ,  —(5) 

and  on  Side  2  it  is  that 

yj +  y2  - 0  •  —(6) 

All  these  relations,  (1)  through  (6),  used  with  the  cavity 
circuit  parameters  shown  in  Figure  3,  will  yield  the  s-plane  resonances 
listed  in  Table  I.  The  numbers  s  =  &¥  jco  correspond  to  resonant  frequencies 
0)  and  to  time  decay  rates  ^"due  to  losses  in  the  cavity  circuit.  Because 
the  circuit  is  very  high  Q  (  100) ,  the  frequencies  u>  are  approximately 

the  same  as  they  would  be  for  a  lossless  cavity  and  lossless  wake. 

The  Table  I.  formulary  can  be  used  to  obtain  expressions  for  the 
circuit  elements  explicitly  in  terms  of  the  resonant  decay  rates  and 
frequencies.  These  are  listed  in  Table  II.  They  are  important,  on  the 
one  hand,  in  relating  measured  circuit  parameters  with  calculated  fre¬ 
quencies  and  decay  rates.  On  the  other  hand,  the  several  Q  formulas  listed 
are  design  criteria  for  the  system  in  that  these  numbers  reflect  the 
manner  of  shifts  of  admittance  zeroes  under  loading  of  the  cavity.  Thus, 
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TABLE  I 

CAVITY  RESONANCES 


Admittance 

Resonances 


s  =  <y+  ju> 


Empty  Cavity 


Low  Zero 

Side  1  Pole 
TMj^q  Terminals 


Side  2  Pole 
TMQ11  Terminals 


High  Zero 


CD  = 
O 


r  C  L 
o  o 


,  The  Measuring  Frequency 


*\/i  1  l 

*1  =  ®1  +  jC°l  =  "  2C~  +  ^  cT  +  lT* 


jj 

o  o  1 


S21  "  ^21  +  ja)21 


-i  .  .  i  w  ~yji  fL  +  L+IL 

=  2Co  (gl  +  q  82}  +  j  *T  Lo  +  Lj_  q  L^ 

s22  =  ®22  +  J*22 

=  -i  ,  .  x  .  sJT  T  i  7T7 

2C  ^81  +  S2^  J  C  +  L,  +  L„^ 


J  AJ  4J-t  Un 

O  O  1  2 


■*  “  <T»  +  :K 


■  ^  [8i  +  (l  +  X]  +V^fc  +  h  +  (1  +  ^] 

“o<a,l<“21>a>22<“3 

0  CO^O^,  <y~22^®3 


INTERRELATIONS 

12212  22  2 
(i  +  “)^22  =  “3  +  and  U  +  =  ^3  + 

(1  +  <5~22  =  ^3  +  \c\  and  (1  +  q)  C* 21  =  <T3  +  *1  0*1 

^22  +  “21  =  “3  +  and  C1  +  O  (“22  "  “21^  =  ^  ^3  "  “l^ 

0*22  21  =  ®"l  and  (  O” 22  ”  21^  =  (  ff" 3  ” 


11 


AEDC-TR-69-32 


TABLE  II 

MEASUREMENT  AND  DESIGN  RELATIONS 


q 


V?2 

klQl 


Measurement  Formulas 


2  2 
to,  -  0)  l 
loo 


0) 


2  2 

ox,  -  Oi  L  ,  L 

-S-2-  f+  <1  +  bf 

CD  L1  q  L2 


2  2 

CD-  -  CD  L 

~~  2"~  <1  +  i>  i: 

CD  M  2 


Q  Formulas 


■2  <yi“o 

. T . - 1  "  ®oLl«l 

CD  -  CD 


<*1 


CD  Q 

O  I 


G  3  1 


CD 


■5iEi  +  (1+?)83 


o’  3  ■  cr  i)  i  i 

oi  q7  (1  +  qJg2 

o  1  ^ 


"2  O'3C0o  ®0L1L2  E2Q2gl  +  (klQl  +  k2Q2)8iI 

k2Q2L2  +  (k1Q1  +  k2Q2)L1 


2  '  2 
®3  ‘  ®1 


-2(  g3  - 

2  2  (A51^j982 

®3  “  ®1 
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for  example,  it  may  be  necessary  to  decide  whether  or  not  a  loaded 

Qi  Q2 

Q  of  the  order  —  or  10  can  yield  readable  measurements  from  the 


system  accurately. 

Figure  4  illustrates  results  given  in  Table  I.  Cavity  and  wake 
are  assumed  lossless  here,  so  that  the  graph  is  one  of  reactance  against 
frequency,  seen  from  either  Side  1  or  Side  2  of  the  cavity.  The  relative 
positions  of  resonances  correspond  to  the  Table  I  listing,  and  0^  is 
assumed  to  represent  either  of  the  admittance  poles  or  ^£2* 

The  frequencies  and  may  be  characterized  very  nearly  as 
modal  resonances.  Strictly  speaking,  pure  TM^q  and  ™oil  modes  dre 
present  only  in  the  empty  cavity,  and  are  only  approximated  to  in  the  filled 


cavity.  In  particular,  the  actual  fields  need  to  have  both  modes  for 
components.  When  the  field  resonates  at  and  ci^ .  the  TMq2q  and  TMqj^ 
modal  components  both  also  resonate,  with  some  coupling  between  them. 

When  the  cavity  field  resonates  at  o>2  (i.e.,  cn2]L  or  o>22)  however,  neither 
modal  component  resonates,  but  each  couples  strongly  to  the  other  mode. 
Thus,  at  o>2^  for  example,  if  the  source  is  at  Side  1,  very  little  energy 
goes  into  y^,  but  most  goes  over  into  y2  and  into  qy^  at  Side  2.  A 
similar  thing  occurs  at  u>22  for  a  source  at  Side  2.  The  frequencies  Ui2 
are  therefore  mode  conversion  resonances.  If  the  TMq2q  terminals  are 
excited,  then  only  a  TMQ11  mode  appears.  If  the  TMg^  terminals  are 
excited  then  only  a  TM^q  mode  appears.  For  an  input  wave  of  one  type, 
a  second  wave  type  is  launched. 

The  measurement  frequency  is  always  below  for  plasmas  and 
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Fig.  4  Loaded  Cavity  Reactance  Curve 
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it  is  always  above  co^  for  non-ionized  dielectrics.  It  could  be  between 

to  and  co,  only  if  a  part  of  the  cavity  held  a  plasma  and  another  part 
1  3 

a 

held  an  ordinary  dielectric.  This,  at  the  same  time,  suggests  a  second 
picture  that  can  view  the  mode  conversion  resonances  cd2.  These  would 
be  resonances  between  the  forward  highly  ionized  wake  and  the  rearward 
much  decayed  wake. 
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SECTION  III 

SPLIT  AND  SHIFTED  RESONANCES  CALCULATION 


Calculation  of  the  resonances  owiyci^  is  done  by  what  can  be 
called  a  frequency  shift  method.  First,  will  be  determined  for 
the  empty  cavity  as  a  simultaneous  resonance  of  its  TM^^  and  tM^q 
modes  in  terms  of  cavity  dimensions.  The  frequency  of  all  actual 
measurement  done  by  the  cavity  is  mQ.  The.  system  will  then  be  perturbed 
by  the  wake,  and  split  and  shifted  resonances  determined  from 

solving  an  electromagnetic  field  boundary  value  problem.  Various 
generalizations  and  details  of  the  solution  —  by  variational  technique  — 
are  discussed  elsewhere  (Refs.  1,  4,  and  5).  Certain  special  adjust¬ 
ments  to  the  technique,  however,  are  argued  here  as  the  calculation  is 
developed. 

The  cavity  is  constructed  as  shown  in  Figure  5(a)  and  is  assumed 
to  carry  only  TM  radially  symmetric  modes.  Circular  cavity  TM  wave¬ 
numbers  are  in‘ general  found  from  the  relation 


£-(rc,n)  2  .  [lit  2  2 

b  d  oo 


—  (7) 


where  b  is  radius  and  d  is  length  for  the  cavity,  where  k2  =  aAi  € 

oo  oHo  o 


is  free  space  wave  number  at  o)q  resonance,  and  1  is  an  integer.  The 

number  p(m,n)  is  an  root  of  an  m-order  Bessel  function,  m  and  n 

integers.  This  applies  to  each  mode  TM 

mnl' 
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+  (!) 2  - k 

b  =  5.520 

1.585 
1=  1.232 


2 

1 


(a) 


(b) 
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Since  circumferential  modes  are  assumed  not  present,  ra  =»  o 
only.  When  operating  TM^q,  then  1  =  o  and  p(0,2)  -  5.520  -  kjb. 

When  operating  TML..,  1  -  1  and  p(0,l)  -  2.405  =  k.b.  The  numbers 
kj^  and  k2  are  wave  numbers  for  the  respective  TM^  and  TMQ1  waveguide 
modes.  Using  this  data  in  relation  (7),  one  obtains  the  simultaneous 
equations 


k 


2 

1 


oo 


and 


—  (8) 


k 


2 

2 


These  can  be  solved  for  b  and  d  as 


—  (8a) 


k  b  =  5.520 

1  OO 


—  (9) 


and 


V 


=  1.110 


—  (9a) 


and,  alternately. 


5  - 


—  (9b) 


2  2  2 

in  which  kj  =  kQo  =  is  used  to  specify  the  frequency  of  operation, 

Since  b . and  d  are  here  real  numbers,  then  the  system  can  operate  in 
these  two  modes  simultaneously. 

The  results  (9) (9a) (9b)  and  equation  (8a)  are  listed  in  Figure 
(5a)  for  reference. 
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Methods  of  solution  of  the  boundary  value  problem  that  would 
characterize  the  responses  of  a  filled  cavity  have  been  discussed 
previously  (Refs,  4,  8,  9,  and  10).  A  variational  approach,  a  Rayleigh- 
Ritz  scheme  such  as  discussed  by  Collin  (Ref.  4),  is  the  one  used  for 
the  .present  cavity  problem.  It  has  the  advantage  that  results  come 
directly  in  terms  of  modes  of  interest  and  couplings  between  them. 

It  is  also  helpful  for  expression  of  results  by  means  of  the  very  small 
changes  in  parameters  (i.e.,  permittivity,  conductance,  resonant  fre¬ 
quency)  that  one  is  concerned  with  here. 

There  are  simplifying  effects  upon  calculations  which  are  due 
to  radial  symmetry  in  the  wake  and  in  the  cavity  modes  that  are  chosen 
for  measurement  but,  most  importantly,  due  to  the  fact  that  the  wake 
causes  only  a  slight  perturbation  of  electrical  properties  of  air. 

More  than  this,  the  approach  which  Collin  shows  is  special,  applying 
just  to  rectangular  guide,  and  would  therefore  need  modification  in 
order  to  try  it  on  circular  guide.  The  modified  and  simplified  equations 
are  derived  below.  The  problem  so  stated  yields  solution  fields  that 
are  approximately  the  correct  ones,  and  leads  to  only  small  errors  in 
the  propagation  constants  and  resonant  frequencies  that  are  to  be  eval¬ 
uated.  Error  incurred  in  derivation  is  discussed  in  Appendix  A. 

Conditions  inside  the  cavity  during  measurement  are  as  shown 
in  Figure  (5b).  The  wake  is  assumed  to  be  cylindrically  shaped  and  to 
have  a  uniform  dielectric  constant  along  its  axial  length  e,  and 

another  uniform  dielectric  constant  M£q  along  its  axial  length  f. 
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The  ratio  e/£  can  be  chosen  arbitrarily,  but  a  special  rule  will  be  shown 
later  to1  specify  e/f  most  favorably  for  accuracy  of  reading  by  the  system. 

The  approach  on  the  problem  is  as  illustrated  in  Figure  6.  First, 
the  two  partially  filled  guides  are  analyzed  separately  by  the  variational 
procedure  mentioned.  They  will  then  be  matched  at  a  common  planar  junction. 
Finally,  resonance  will  be  said  to  occur  when  nodes  (on  either  side  of 
the  junction)  are  separated  by  total  length  d.  Frequencies  found  this 
way  are  the  modal  resonances  and  (see  Figure  4.)  The  mode  conver¬ 
sion  frequencies  on  the  other  hand,  are  not  resonances  in  the  same 
sense  of  the  procedure  used  to  find  and  cd^.  However,  these  can  be 
found  easily  enough  from  circuit  interrelations  listed  in  Table  I. 

Now  let  K(r)  be  the  radial  function  giving  relative  dielectric 


constant,  i.e. 

K(r)  =  K  ,  o^r<a 

- (10) 

=  1  ,  a<r  ^  b 

Then  for  this  cylindrical  region  the  pertinent  Maxwell's 

be  written 

equations  can 

V  x  |H  =  j“>€GK(r)  |E 

1 

1 

1 

t-* 

l-» 

Vx  lE  = 

- (lla) 

V*  |H  =  0 

- (lib) 

V*  K(r)  |E  =  0 

- (11c) 

20 


Fig.  6  Steps  of  the  Resonance  Calculation 
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In  order  not  to  labor  familiar  parts  of  the  derivation  process,  it 
is  noted  first  that  only  TM  waves  are  of  interest,  and  second  that 

'ft 

|E  and  |H  can  be  given  in  terms  of  potentials  77Tand  0  by 

|H  =  &€0SJx'Tir  - (12) 

and,  where  k2  =  0)2m  €  , 

0^0  0 

2 

l E  =  ko7F+y  0  —(12a) 

But  for  a  factor  K(r),77T  would  be  an  electric  Hertzian  potential.  These 
two  relations  can  be  inserted  into  the  first  listed  Maxwell's  equation 
(11)  to  obtain  the  ordinarily  crucial  statement 

y27T+  k2K(r)7F=  V(V  •W)-K(r)^  0  —(13) 

A  solution  of  this  is  what  is  needed.  Auxiliary  relations  will  now  be 
discussed  which  can  simplify  the  steps  of  the  calculation. 

Suppose,  on  the  right  hand  side  of  equation  (13),  that 

V<  V'TTT)  -  -  0  —(14) 

This  condition  is  usually  made  when  specifying  potential  in  homogeneously 
filled  waveguide.  The  somewhat  changed  consequences  of  using  (14)  for 
the  present  partially  filled  guide  need  to  be  characterized.  Putting  (14) 
with  (13),  then  /ff  found  from 

V2? 7T+  k2K(r)7W=  -  ^(r)  -  lj  y(  y  )  —(15) 


See  Collin  (Ref.  4)  or  Stratton  (Ref.  ll)  for  more  thorough  treatment  of 
potentials . 
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would  be  the  field  required  to  solve  the  problem.  In  this  case,  the 
associated  electric  field  becomes 

ie  =  k^-nr+V  <V  ,7F>  *  —(i6) 

having  used  equation  (14)  to  substitute  for y  0  in  equation  (12a). 

Continue  this  supposition  by  letting77Ttake  the  form  of  a  z-directed 
E-wave  potential 

77 T  =  ^  exp^z  -“"(17) 

where3jf  =  3jf(r)  fits  zero  boundary  conditions  on  the  guide  walls  of 
Figure  6,  and  where  the  propagation  constant^  is  a  number  still  to  be 
determined.  Used  in  equations  (12)  and  (16),  this  gives  for  the  electro¬ 
magnetic  fields  the  statements 

|H  =  -jas  exp^z  —(18) 

and 

/E  -  f(k*  +y2)^^+yr  |^j  exPyZ  ,  —(18a) 

where  £  are  radial,  circumferential,  and  axial  respective  unit 

vectors.  Plainly,  with  and  continuous,  these  |E  and  |H  fields 

fit  correct  transition  conditions  across  the  boundary  between  the  cylindri¬ 
cal  wake  and  air  inside  the  guide  in  Figure  6,  as  well  also  across  the 
planar  junction  between  the  two  wake  sections  K  and  M. 

This  same  potential77ysubstituted  into  equation  (15)  gives  the 
two  component  vector  statement 
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[7  77  (r  lr) +  +72)K(r)^j  S' 

+  |K(r)  -  IJy^r  =  0  ,  —(19) 

which  is  incorrect,  but  which  is  small  in  error  even  assuming  the  £ 
vector  component  zero.  The  purpose  of  doing  so  is  that  solution  of 
(19)  becomes  simplified  thereby  and,  moreover,  that  electric  field 
expression  in  terms  of  the  potential  function  thereby  also  simplifies. 
Thus  the  equation  that  needs  to  be  solved  here  is 


7  57  (r  ir)  +  K<r>  K  +T2>  0  .  —(20) 

and  the  expression  of  |E  and  |H  fields  is  done  by  formulas  similar  to 
those  used  with  homogeneously  filled  guides.  The  character  and  size 
of  error  incurred  in  forgetting  the  r-component  of7JT  by  using  (14)  is 
discussed  in  Appendix  A. 

A  suitable  form  for"^,  for  solution  to  (20) ,  can  be  as  a  sum  of 
purely  radial  mode  functions  of  the  circular  waveguide;  i.e-. , 

y-X  %  'f'n  ■  —(21) 

n 


where  the  a^  are  amplitude  coefficients  and  where  the  are  the  unit 
functions 


1  “\/r j  (k  r) 

■jt  =  °  n 

Tn  b  J^b) 


—  (22) 


Here,  Jq  and  are  zero  and  first  order  Bessel  functions.  The  wavenumbers 
kn  are  found,  as  specified  in  equation  (7),  such  that  knb -is  an  n—  root 
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of  the  zero  order  Bessel  function.  The  property  that  makes  \^/n  a 
unit  function  is  that 

ll^nll  2  =  J  ^/rdr  =  1 

Associated  with  equation  (20)  is  the  quadratic  relation  (Refs. 

4  and  5) 

J\zr\  rdr  -  J  K(r)[f0l02  +  ko]  Jvjrfrdr  =  0  ,  —(23) 

2  /»  y  2 

where  the  function  rojo*  has  the  stationary  value J  whenever  the 

correct  potential's  used.  By  being  stationary,  percent  error  in  cal- 
2  2 

culating'y  +  kQ  this  way  is  at  worst  proportional  to  percent  mean 

squared  error  in  This  is  error  due  to  approximation  of  fields  by 

a  finite  number  of  modes.  It  is  not  the  same  as  the  operator  approxi- 

* 

mation  for  (20)  mentioned  earlier.  In  fact,  it  is  some  orders  less  an 
error  than  that  in  (20) . 

An  alternate,  more  convenient,  form  for  the  above  relation  (23) 
can  be  given.  Using  the  fact  that  K(r)  has  constant  values  inside  and 
outside  the  wake,  one  can  write 

l  |srf  '<*  -  [roio2  +  *3  [I*1  +  CK-D  /.|tf  «*]  -  o.  — 

where  b  is  radius  for  the  guide  but  a  is  radius  for  the  wake.  Assuming 
to  have  the  linear  form  shown  by  relation  (21),  then  integral  terms 
in  (23a)  correspond  to  matrix  operations  (see  Collin,  Ref.  4)  listed  in 

*This  can  be  laid  to  the  fact  that^is  twice  differentiable,  having 
higher  modal  components  which  are  very  slight. 
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TABLE  III 

OPERATOR  AND  QUADRATIC  FORM  CORRESPONDENCES 


1 M2'*' 


T.  w 


I  \m  -  <*• 


Lm 


rdr  <3~ 


where  A.  .  =  A. .  =  /  W/.  w/  rdi 

1J  J1  Ja  ^  ^ 


■  ^  a/Tw> 

and  n  -  b  j^K.b) 


o.. 

k? 


r*ll 

*12 

*13 

*21 

*22 

*23 

*31 

*32 

‘  *33  . 
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Table  III. 

Using  these  results,  and  assuming  only  the  two  effective  modes 
TMq2  (for  n  =  1)  and  TM^  (for  n  =  2) ,  then  solution  of  the  variational 
problem  is  equivalent  to  solution  of  the  matrix  equation 

(*i  -  O'2  +  k*>  t1  +  «-i>^  -<72  +  k2)«-i)^2  ' 

V  -  <72  +  k2)(K-l)A12  ^  -  <y2  +  k 2)  [l  +  (K-l)^ 

—  (24) 

The  determinant  of  the  matrix,  set  to  zero  and  solved,  yields 


In  (25),  were  the  propagation  constants  %(-  0)  and^  (*  J  ^) known, 

then  k  =  k  ~  and  k  =  k  ,  respectively  would  become  known,  and  then 
o  o3  o  ol 

2  -12  2  -12  ft 

cd3  =  (no  €  )  kQ3  and  =  (|aQ  €  Q)  kQ1  respectively  could  be  deter¬ 

mined.  The  resonances  must  next  be  found  by  marching  this  guide  (wake  K) 
with  the  guide  having  wake  M,  and  then  closing  off  the  ends  so  that  it 

•jc 

It  turns  out,  as  the  1-2  vs  3-1  subscripts  indicate,  that  the  it  resonance 
(TMQ1 ..)  has  a  smaller  frequency  change  than  the  0  resonance  (TM^g)  .  This 
is  intuitively  sensible,  knowing  that  the  waveguide  is  a  high-pass  filter 
whose  cutoff  frequency  is  raised  by  the  wake.  It  means  group  velocity  is 
higher  for  the  cutoff  TM«2  wave,  and  this  means  a  higher  frequency  change 
to  reestablish  resonance  for  the  TM^q  mode  than  for  the  TMqj^  mode. 
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resonates  only  over  the  length  d,  as  suggested  in  Figure  6. 

With  very  small  error  in  approximation  (even  for  a  wake  Q  of 
2 

100)  ,  th expression  (25)  can  be  put  more  conveniently  by  formulas 
listed  in  Table  XV.  The  same  adjustment  worked  for  the  second  waveguide 
(substituting  M  for  K)  yields,  with  the  same  qualifications,  results 

2  2  jr  2 

for  h^  =  0  and  h2  =  -  ^  also  shown  in  Table  IV. 

Table  V  shows  potential  field  functions ^ that  can  exist  in  the 
two  waveguides,  in  terms  of  the  basic  radial  mode  functions  [J/  .  It  is 
linear  combinations  of  these  fields  that  one  would  generally  use,  depend¬ 
ing  on  particular  guide  end  and  guide  junction  conditions.  The  coeffi¬ 
cient  matrices  listed,  because  of  the  interrelations  shown,  are  not 
orthogonal.  However,  because  of  small  size  in  the  cross  diagonal  elements, 
they  are  approximately  orthogonal.  These  are  the  same  coefficients 
needed  to  solve  equation  (24).  The  designations  TMq2  and  TM^  are  also 
only  approximate  since  now  both  modes  are  present  in  each  field. 

Match  at  the  junction  of  M  and  K,  indicated  in  Figure  6,  can  be 
made  in. terms  of  transverse  electric  and  magnetic  field  components  derived 
from  the  ^potentials .  Thus,  from  equations  (18)  and  (18a)  and  from 
Table  V, 

=  ^sh  7l(z  +  e)  +  A,  5^  Coshy2(z  +  e)  —(26) 
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TABLE  IV 

LOADED  GUIDE  PROPAGATION  FUNCTIONS 


™02 


Relation  to  End 
Results 


2  2 

y.  +  k  2 

1  2  °  =  1  -  (K-l)&u  +  ^K_1> 


r  2  /kid\2  2 " 

Ki  +prj  ^2 


7i  * 0 


,2  ,2  2  * 
k  =  kn-  =  03,u  ^ 
o  03  3  o  o 


™01 


2  2 

y2  +ko 


=  i 


('■DV  +  «-l)2[^2  '('ll  ^ 


™02 


?‘-0  2 

ko  *  k0l"  4o€, 


Relation  to  End 
Results 


2  2 

h,  +  k 


1  '-IT  •  i  '  <M'1)Au  +  (lf'1)2[4  +(-^  4] 


bii0 

2  2 
k  =  k 
O  03 


Voe, 


TM, 


01 


2  2 
h2+ko 


=  1 


“  +  (M'1)2p22"  (“rj  *12 


2  .  a2 
2  "  “  6) 

,2  ,2  2 

k  =  k.  =  03-li  € 

o  01  l^o  ^o 
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TABLE  V 

LOADED  GUIDE  POTENTIAL  FIELDS 

’/'i-'k  FROM  TABLE  III ;  y,  h  FROM  TABLE  IV. 


^1K  ^*11^1  +  a2lVV 


^1M  "  *bll^l  +  b21^2^ 


*2k  =  <a12^1  +  a22^2) 


^2M  =  (b12^1  +  b22^V 


811  ®I2 


La21  a22 


bll  b12 


b2I  b22 


a21  /kld  2 


b21  -  kld  2 

b^  =  —  (M-l)A^2 


S‘ -(¥)■«• 


,2  ,  2 

kla12  "k2a21 


k2d\2 


a  I  (M“1)A^2 


klb12 


-k2b21 


®11  ,'®22 


bll  =  b22 
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and 


tK 


7iAiS?ksi“h'X(2 


^')Lr0v 

+  •)  +y2A2  sr3  sinhy2(i 


+  e)  —(26a) 


give  transverse  fields  in  the  K  guide,  while 

is?;  \  -  h  ^ Cosh  hi<*  -  «  +  b2  sr^ Cosh  V*  -  £>  — (27> 


and 


EtM  “  hlBl  Sr 


dl?lM 


dy. 


Sinh  h^z  -  f)  +  h.^ 


2M 


Sinh  h2(z 


f) 


—  (27a) 


give  transverse  fields  in  the  M  guide.  At  the  junction,  in  these 
expressions,  one  sets  z  =  0,  =  0,  and  E^  -  E^  =  0.  In  addition, 

Z\b 1  Zdj 

since  r1 —  and  r-1 —  are  still  orthogonal  (and  taking  orthogonal  components) , 
or  or 

one  is  led  to  the  four  equivalent  statements 


/  ancoshyie 

a12Cosh  ye 

-b^Cosh  h^f 

-b^Cosh  h2f  \ 

/M 

a21Cosh  ^e 

a22Coshy2e 

-b2iC°sh  h^f 

-b22Cosh  h2f 

A2 

l  7iansinh7ie 

y2ai2sinh72e 

h^b^Sinh  h^f 

h2b12Sinh  h2f 

K 

Vja21Sinh71e 

inline 

hlb21Sinh  hl£ 

h2b22Sinh  h2f' 

\»  J 

put  in  matrix  form. 


—  (28) 
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The  determinant  of  the  (28)  matrix,  set  to  zero,  yields  the 
characteristic  equation  for  the  junction 


h,h„  tanh  h,  f  tanh  h„f 


12 


yiYl  t®nh'yie  tanh-y2 


h  tanh  h..  f  h  tanh  h„f 
- —  +  1  +  ~  - r--—  '  4-  - ~ — 


tanh^e  -y2  tanh'y2e 


+  (CUC22  -  1) 


•  V*.!  % 


[h^  tanh  hjf  h2  tanh  l^fl 
yl  tanh  +  y[  tanhy2ej 

-] 


tanh  h..  f  h„  tanh  h_f 

1  +  ~  r— ,  ■  2  1  =  0 


17*  tanh/y2e  ’  y^  tanh'y^e 

where  the  C  multipliers  are  found  as  coefficients  of  the  matrix 


—  (29) 


By  using  (30)  and  results  in  Tables  IV  and  V,  one  can  show  the  last  two 
bracketed  terms  in  (29)  to  have  value  proportional  to  (C^^p  •  This 
is  by  some  orders  smaller  than  previously  mentioned  errors.  Ignoring 
these  terms,  then  the  junction  characteristic  equation  (29)  can  be  written 
in  the  much  simplified  form 


[h^  tanh  h^f  ”1  tanh  h2f  "I 

Y  ta"hy1e  +  XJ  |yl  “nhy2e  +  j ' 0 

If  one  takes  advantage  of  the  fact  that  propagation  constants  of  the 
loaded  guide  are  but  slightly  changed  from  those  of  the  enjpty  guide. 


32 


AEDC-TR-69-32 


then  (31)  can  be  rewritten 


[*  ♦  *  0  *]  *-(VMr]a}-° 


-(31a) 


These  last  two  expressions  of  the  junction  characteristic 
equation  point  to  an  interesting  property  of  the  cavity:  namely,  that 
the*y^  and  h^  modes  affect  only  each  other,  as  also  do  the  an<*  hj 
modes.  The  chief  perturbing  influence  on  behaviour  of  the  cavity,  in 
other  words,  is  the  wake  ionization  in  front  relative  to  ionization 
behind.  Qualitatively,  this  is  good.  It  suggests  it  is  only  a  mild 
assumption  —  even  at  high  front  to  back  ionization  ratios  --  to  hold 
that  the  wake  can  be  characterized  by  a  sudden  change  in  its  dielectric 
constant.  It  is  not  to  say  there  should  not  be  such  coupling,  but  without’ 
coupling  there  is  less  that  needs  explaining.  The  matter  is  argued 
briefly  as  follows. 

Use  of  equation  (31)  is  equivalent  to  the  assumptions 

a12  ‘  a21  -  b12  '  b21  '  °  a“d  aU  =  a22  -  bll  ‘  b22  "  1 
in  the  statement  of  equation  (28).  In  this  case,  that  homogeneous  ex¬ 
pression  can  be  rewritten  as 


/ 


Cosh'^e 
0 

7]S  inline 

\° 


Cosh^e 


'Vsinh')*' e 


-Cosh  h^f 


h^Sinh  h^f 


-Cosh  h^f 


h^Sinh  1*2  f 


=  0 
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It  is  plain  that  solution  vectors  for  (32)  must  look  either 


—  (33) 


There  are  no  other  solution  vectors  for  (32) .  These  forms  (33)  agree 

with  the  qualitative  suggestion  of  cavity  behaviour  just  made. 

The  junction  characteristic  equation  (31a),,  substituting  ~y and 

h  expressions  from  Table  IV,  can  be  solved  for  kQ1  and  kQ3.  For 

^h^,  one  sets  kQ  =  kg3  to  solve.  For^h^,  one  sets  kQ  =  k^  to  solve. 

Results  (with  negligible  additional  error)  are  listed  in  Table  VI  in 

2  2  2  2  2 

several  forms.  The  relations'  k  =  k  =  o>  u  6  ,  k„„  =  cd,u  (  ,  and 

1  oo  o  o  o’  01  lHo^o’ 

2  2 

k03=  0>3|jio €  were  also  used  for  the  expressions  in  Table  VI.  It  is 
interesting  to  note  here  that  the  number 


.K-l  =  |  (K-l)  +  |  (M-l)  —(34) 

is  mean  susceptibility  while  the  number 

1  (K-l)2  =  |  (K-l)2  +  |  (M-l)2  —(35) 

is  only  a  formal  mean  square  susceptibility.  That  is,  in  the  first  place, 
(35)  gives  information  derived  solely  from  the  mean  value.  In  the  second 
place,  both  (34)  and  (35)  retain  complex  number  information  about  wake 
susceptibility.  Thus,  while  the  single  mode  cavity  is  well  known  to  give 
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TABLE  VI 

CAVITY  WAVENUMBERS  AND  RESONANCES 


2  2.2  .2 
"l  -  “o  k01  -  kl 


2 

01 

0 

'  k! 

—  ss 

2 

“3 

2 

-  03 

O 

k03  ‘ 

2 

03 

0 

'  k? 

2 

!2l 

2 

-oil 

k03  ' 

k01 

*  *22  +  W’1*2  [^)  [4  -(¥)  4]i2H 


1  -  *  (K-1’*!!  +  «-»2  [4  +(¥)  4] 


=  -(FT)  fii  *  (^) 2  *22] 


High 


4  .  ,4 


2  T  2  /^2\  ^  2  ^1  +  ^2  /^2^\  2  2  ”1 

<*-»  [firp]  4  +  -^  ft]  4J 


kk_Ji  .  .  (^2  .  (K-2)2  ^-f^‘4] 


=  “  (K'1)^2  +  <K-X>  P22  -  - 


K-l  =  t  (K-l)  +  4  (M-l)  Mean  Susceptibility 

a  a 


(K-l)2  =  4  (K-l)2  +  4  (M-l)2  Formal  Mean  Square  Susceptibility 

u  O 
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an  averaged  result ,  the  two  moded  cavity  is  here  seen  to  give  both  an 
averaged  and  a  kind  of  averaged  square  result.  This  property,  it  will 
be  shown  later,  provides  help  toward  improvement  of  accuracy  of  system 
measurements  and  toward  an  increase  in  information  available  from  the 
system. , 

There  is  error  in  the  assumption,  for  results  in  Table  VI,  that 
K  and  M  can  be  complex  numbers.  However,  it  is  small  error  —  of  the 
same  type  and  amount  as  that  in  the  approximation' of  equation  (19)  by 
equation  (20)  —  and  therefore  accounted  negligible.  When  K  and  M  are 
complex,  then  the  frequencies  a^,  a^,  and  cd3  are  also  complex.  Indeed, 
say  such  that 


s  =  C7+  jcu 


—  (36) 


indicates  the  replacement  necessary  for  these  respective  frequencies. 

Now  consider  the  effect  on  Table  VI  formulas.  One  has  it,  for 
example  when  jo^ - £>s1  =  +  jo^  and  jov — J>sq  =  jcDQ,  that 


i.  -  s  s.  +  s 

1 _ O  1  o 


o  o 

<r 1  +  31  +  8< 

jet  s 


<Ti  si 


S.  +  S  CD  -  £15  S  '+  S 


JCD  S 
J  o  o 


°+  1 


o  1 


CO 
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<Fl 


2  + 


0),  “  0)  CO,  +  to 

1 _ o  1  o 


CO 

o 


CO 

o 


CO,  -  CO 

1  o 


CO 


•  Cl 

y°n 


—  (37) 


o 

The  approximation  is  a  good  one,  since  it  falls  upon  frequency  differences 
and  not  upon  absolute  frequencies.  The  same  treatment  shows  also  that 


o 


and  that 


2  2 
s3  "  S1 


2  2 
co3  -cox 


CO 


2 (cr3  -crj 

3“ 


—  (39) 


By  comparing  these  frequency  and  decay  rate  expressions  (37)  (38)  (39) 
with  equivalent  ones  in  Table  II  and  in  Table  VI,  one  sees  there  are 
formulas  now  to  relate  cavity  equivalent  circuit  measurements  with  cavity 
electromagnetic  field  calculations.  These  are  listed  in  Table  VII.  One 
can  notice  these  expressions  are  linear  in  the  susceptibility  averages, 
and  so  solve  for  them  in  terms  of  known  measurable  and  calculable  quanti¬ 
ties  of  the  system,  as  also  listed.  The  last  formula  listed  is  a  design 
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TABLE  VII 

CIRCUIT  MEASUREMENTS  RELATED  TO  MODAL  CALCULATIONS 


Continued  next  page 
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TABLE  VII  (Concluded) 


|M2  fn  .  (Wl2 

'N  ^ 


k,b  =  5.5201,  1(0]  =  0.3403,  kjb  =  2.4048,  J,(k-b)  =  -0.5191 
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figure  of  merit  for  the  system.  It  is  derived  assuming  K-l  small  and 


(K-l)  negligible.  What  it  asks  is  that  the  design  be  done  so  that  the 
two  measured  admittances  be  roughly  the  same  order  of  size.  Failing 
that,  one  could  need  to  inquire  into  the  accuracy  of  measurement  of 
them.  Table  VII  gives  the  final  results  of  this  section.  There  remains 
just  the  special  choice  of  e/f  cavity  axial  division  to  give  and  to 
interpret. 
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SECTION  IV 

OPTIMIZED  CAVITY  MEASUREMENT 

The  formal  averages  in  Table  VII,  once  a  ratio  is  specified, 

can  let  one  calculate  an  M  in  the  f  section  and  a  K  in  the  e  section 

of  the  cavity  (i.e.,  by  solving  equations  (34)  and  (35))-  However, 

every  choice  ^  lets  some  corresponding  K  and  M  be  calculated.  One 

could  ask  whether  there  is  not  some  optimum  choice  which  should  be 

made  for  a  preferred  K  and  M  pair;  in  particular,  to  be  done  without 

requiring  additional  measured  data.  A  preferred  -  ratio  can  be  that 

which  locates  the  junction  where  actual  electric  susceptibility  is  the 

same  as  measured  average  K^.  This  is  all  the  more  so  a  good  choice 

because  it  also  means  that  K-l  and  M-l  are  actual  averages  over  the 

0 

respective  e  and  f  cavity  sections.  Most  importantly,  this  —  marks 
the  wake  at  a  particular  distance  behind  the  model  such  that  subsequent 
graphs  are  then  precisely  tied  to  model  location. 

A  positive  choice  can  be  made.  The  question  raised  is  perhaps 
not  such  a  simple  one  as  the  approach  in  answer  to  it  reflects. 
Nonetheless,  it  is  a  step  in  the  right  direction.  Indeed  the  suggestion 
that  is  made  is  such  as  should  be  applied  to  all  data.  The  matter  is  a 
physical  one.  Essentially,  it  attempts  to  make  the  junction  (the 
sudden  susceptibility  change)  appear  as  unobtrusive  as  possible.  In 
doing  so,  it  matches  the  physical  character  of  the  wake  at  the  same 
time.  One  result  is  to  allow  higher  M-l  to  K-l  ratios  to  be  read  than 
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for  the  case  e  =  f  =  y,  but  this  is  a  property  of  any  move  away  from 
a  half-half  split  in  cavity  axial  length. 

The  matter  is  stated  as  follows.  Consider  -two  slices  of  the 
wake  as  shown  in  Figure  7  (or  even  smaller  constant  volume  elements). 
The  energy  difference  between  the  Ae  and  Af  sections  can  be  much  larger 
at  the  junction  than  if  the  slices  were  moved  together  anywhere  else 
along  the  wake.  By  comparison,  if  this  were  the  real  wake  without 
junctions,  differential  energy  between  the  slices  would  still  vary  as 
the  pair  were  moved  abcut,  but  variation  would  be  by  a  second  order 
of  slice  size  —  not  by  a  first  order.  One  can  say  that  large  energy 
differential  at  the  junction  corresponds  to  first  order  error  which 
should  be  reduced. 

The  energy  in  the  wake  is  almost  entirely  electric.  In  the 
slices,  it  is  proportional  tc  (K-l)Ae  in  the  e  section,  and  to  (M-l)Af 
in  the  f  section.  Proportionality  is  the  same  for  both  sides  since 
both  sides  have  approximately  the  same  field  strengths.  Therefore, 
the  difference 


(K-l)Ae  -  (M-l)Af  =  A Q,  —(40) 

is  the  thing  to  minimize.  Now  equations  (34)  and  (35)  can  be  solved  to 
give 


K-l 


—  (41) 
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Fig.  7  Volume  Element  in  Differential  Energy  Test 
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and 


M-l  =  K-l  -T'vf-  "V(K-l)2  -  (K-l) 2 


—  (42) 


Substituting  into  relation  (40)  gives 


a£=  K-l  Ad 


—  (43) 


where  Ae  +  Af  =  Ad  is  constant.  Differentiating  A^with  respect  to 
e 

—  gives  the  minimizing  condition 


Ae  e 
Af  "  f 


—  (44) 


This  result  put  into  (40)  yields  finally  the  junction  location  con¬ 
dition 


(K-l) e  =  (M-l) f 


—  (45) 


The  energy  that  would  pour  through  the  junction  because  of  sudden 
susceptibility  change  is  thereby  reduced  at  this  position.  Looked  at 
from  another  point  of  view,  condition  (45)  holds  that  there  is  the  same 
weight  of  ionized  material  in  the  forward  part  as  in  the  rearward  part  of 
the  cavity. 

The  equations  (45) (42)  and  (41)  (or  (45) (34)  and  (35)  )  can  be 
solved  to  put  e  and  f  explicitly  in  terms  of  measurements,  and  are  all 
collected  and  placed  in  Table  VIII.  These  results  constitute  the  end 
point  of  the  measurement  procedure  discussed  heretofore. 


44 


( 


K-l  =  K-l 
M-l  =  K-T 


e  _  M^l 
f  =  K-l 

f  =  K~2_ 
e  M-l 

where  A 


TABLE  VIII 

OPTIMIZED  MEASUREMENT  RELATIONS 


V 1-  A2 


A/l  -  a2 


\  Junction  Location  Condition 

Is  ->/7TT] 2 

I 

....  2  - 

g-2  _  (K-l?  ..  and  Where  (X  -  T) and  (K-l)  are  found  from 

(K-l)2 


measurements  given  in  Table  VII 
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APPENDIX  A 

ERROR  IN  APPROXIMATING  THE  WAVE  EQUATION 


The  manner  by  which  error  is  considered  here  rests  on  the 
facts  that  the  loaded  guide  is  but  a  slightly  perturbed  empty  guide, 
and  that, all  field  components  are  but  slightly  perturbed  by  about 
the  same  proportional  amount  from  empty  guide  fields.  In  this  view 
more  particularly,  the  correction  field  is  not  just  small  but  is  at 
the  same  time  smooth.  The  modal  amplitude  coefficients  decrease  rapidly. 
Therefore,  even  a  field  approximated  by  one  or  two  modes  can  serve  for 
a  test  of  the  wave  operator  in  (15)without  causing  serious  effects  due 
to  the  finite  mode  approximation.  It  can  be  assumed  because  of  this 
that  loaded  guide  operators  work  on  fields  very  nearly  the  same  as 
those  of  the  empty  guide,  and  yield  valid  estimates  for  error  due  to 
operator  approximations. 

The  following  property  of  stationary  eigenvalues  \  of  a  symmetric 
operator  A  is  useful.  The  error  aX  is  of  the  second  order  and  can  be 
written 


<y.  a  y)  -\||yn  2 

llx  +  y  II 2 


—  (46) 


where  x' is  an  exact  field,  x  +  y  an  approximate  field,  and  y  the  error 
component.  If  A  has  an  eigenvalue  X'  very  close  to  X.  and  if  y  can  be 
held  to  the  corresponding  eigenvector  for  X 1 »  then  the  error  (46)  can  be 
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rewritten  as 

*X-<X'-X>  , 113,11  ~2  •  -(47> 

thereby  bounding  corrections  put  on  In  particular,  this  is  good 
because  only  some  lowest  eigenvalues  of  A  have  interest  here,  and  these 
are  well  enough  separated  between  thems elves  and  from  higher  characteristic 
values . 

Now  the  potent ialTTT that  solves  equation  (15)  can  more  generally 
have  the  form 


w  -f  £  +  *  y]  «p(  yo  .  ---<«> 

where  <^and  V  are  approximately  proportional .  In  the  empty  guide  as 
in  the  completely  filled  guide, Cj)  and ^ can  be  identical.  If  (48)  is 
inserted  into  the  vector  equation  (15),  then  the  pair  of  scalar  equations 


1  d_ 

r  dr 


+  K(r)  (y2  +  k^)Tjr  +  (X(r) 


and 


«<r,  -  *  «*>  a  -  o  ~<«o 

can  be  put  down.  If  it  is  taken  that”Vjf  =  and  y'P  •  where  \fj  is  a 
unit  Bessel  function  from  Table  III  and  where  corrections  to ^ and  (J) are 
small  and  nearly  proportional,  then  the  equations  (49)  and  (49a)  can  be 
put  in  the  equivalent  matrix  form 


49 


AEDC-TR-69-32 


(-k2  +  (72  +  k2)  [I  +  (K-l)Aj 
(K-l)7kA 


-(K-l^kA 

(-k2  +  k2)  (l  +  (K-l)A]+7^l  yi 


I-  0. 


—  (50) 


Here,  k  and  A  are  kj.  and  An  for  the  TM^  mode  and  k£  and  A^  for  the 
TMQ1  mode.  It  will  not  be  necessary  to  specify  k  and  A  to  correspond  to 
either  one  or  the  other  mode.  The  same  treatment  would  serve  for  each  of 
them. 


The  determinant  for  (50),  set  to  zero  and  solved,  yields  the  two 
propagation  constants 

7«l>)2  =  -k2  +  k2  .  k2  (K-l)A  —  (5i; 

and 


yOP)2  =  "ko  +  »  —(51a) 

associated  with  the  <j>and 3J/ potentials  respectively;  both  very  nearly  zero 
This  is  in  contrast  to  what  happens  with  loadings  for  travelling  wave 
tubes.  There,  one  needs  to  separate  two  modes  30  as  to  be  able  to  use 
the  slower  wave  (Chodorow  and  Susskind,  Ref.  3),  Here,  two  modes  are 
needed  in  order  to  determine  a  maximum  phase  change  for  the  system  under 
loading.  The  phase  change  can  be  determined  from 

y2  =yojo2  +y«i»2 

=  -k2  +  k2  +  k2(l-K)A  ,  ...(52) 

for  which  the  following  three  comments  apply.  Firstly,  both<J> andTjr 
potentials  contribute  to  the  split  and  shift  of  resonance  discussed 
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in  Sections  II  and  III.  Secondly,  the  result  (52)  is  derived  by  working 
through  the  junction  matching  procedure  given  in  Section  III.  Thirdly, 
the  approximation  (15)  holds  even  giving  more  modal  components  to  the 
potentials  <j>  and  The  result  (52)  is  thus  taken  for  a  lower  bound 

on  values  ^ . 

For "y($)2  and  “/( JO2  in  (51)(51a),  one  finds  from  (50)  either 

|2|  or  ||  =  VSToI  ,  —(53) 

depending  on  the  potential  mode  figured.  This  means  that! 7T ,  in  terms 
of  its  component  amplitudes  can  look  like  either  the  vector 

£  ~\/( K-1)A  +  £  (54) 

for'yC'JO  >  or  else  tlie  vector 

r  +  z-\/CK-l)A  (54a) 

for  7«I» .  The  vectors  (54)  and  (54a)  would  have  the  same  magnitude 
error,  were  the  smaller  component  overlooked  in  each  of  them*  However, 

2 

it  is  the potential  which  is  the  more  useful  one.  The  reason  for  ^ 0|0 
so  small  in  (51a)  is  due  to  coupling  between  ¥  and  (£)in  (50)  and  in 
(49) (49a) .  Without  coupling,  the  solution  to  (49)  will  yield  a  propagation 
constant  bounded  by  (52)  and  by  an  upper  bound  to  be  discussed  now. 
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Equations  (49) (49a) ,  without  arguing  this,  can  be  put  suitably 
to  test  the  vector2(54)2in  the  error  criterion  (47).  Doing  so,  letting 

-  1  and  a\  the  error  in  this  quantity,  then 


be  a  value  of 


u:  +y 


A.  IA  A.li  +  |(K-i)a  | 

=  | X'  -\||(k-i)a| 


- (55) 

requires  just  an  estimate  of  the  eigenvalue  separation  |X'  -  Xi*  A 

lower  bound,  already  given  by  (52),  an  upper  bound  is  needed. 

-nu.  y2  +  k : 

This  is  shown  from  consideration  of  - - = -  -1  found  from  three 

k 

possible  wave  operator1  approximations .  The  three, 

Vt  2  +  K(r)(y2  +  k2)  ,  - (56) 

Vt  2  +y2  +  K(r)k2  ,  —(56a) 

K(r)yt2  +y2  -H  K(r)k2  ,  —(56b) 

all  working  on  radial  functions,  are  not  hard  to  interpret.  All  are 

correct  for  an  empty  guide.  Only  (56a)  is  correct  for  a  filled  guide.  A 

comparison  of  their  effects  on  \  is  shown  in  Figure  8,  where  the  scale 

2  2 

(K-l)A  is  proportional  to  only  when  is  small. 

The  operator  (56b)  corresponding  to  the  uncoupled  ^potential,  yields 
X always  too  small.  The  operator  (56),  corresponding  to  the  uncoupled*^ 


* 

These  are  rearrangements  using  operator  formalisms. 
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potential,  provides  a  worst  case  upper  bound  on  X  at  the  measurement 

w  ■  * 

frequency  to  .  The  operator  (56a)  provides  a  next  worst  case  upper 
°  2 

bound  on  X  for  small  at  the  frequency  of  loaded  cavity  resonance. 

The  best  case  curves,  also  at  loaded  cavity  resonance,  are  not  used 
for  the  bound  evaluation  but  are  meant  to  illustrate  the  effect  on\ 
of  coupling  the  ^>and  ¥  potentials  as  in  equations  (49)  (49a)  ....'From 
equation  (52),  and  from  the  best  case  curves,  X would  be  the  sum  of 
the  and  curves . 

From  the  worst  case  curves,  the  upper  bound  for  \  yields 


"'32 

(K-1)V_... 

K 


—  (57) 


which  can  be  suitable  for  larger  wake  diameters 

-  *  **' 

case.  For  the  assumed  present  wake  size  (small 
worst  case  curves,  the  upper  bound  for  Xy*-elds 


2 


than  in  the  present 

2 

) ,  from- the  next 


a\  =  |  (K-l)A  |  3 


- (58) 


The  error  estimate  (58)  implies,  if  equation  (15)  is  used  to  represent 
the  loaded  cavity,  and  if  a  wake  with  Q  (=  |K  -  1|  “1A“1)  about  100  is 
measured  for  example  (enough  to  short  the  system  out),  then  calculation 
error  could  limit  the  determination  of  wake  average  location  to  an  interval 
of  0.1  inch  out  of  a  10  inch  long  cavity. 
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